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Abstract 

Leaf equivalent water thickness (LEWT) is an important indicator of crop water status. Effectively monitoring the water 
status of wheat under different nitrogen treatments is important for effective water management in precision agriculture. 
Trends in the variation of LEWT in wheat plants during plant growth were analyzed based on field experiments in which 
wheat plants under various water and nitrogen treatments in two consecutive growing seasons. Two-band spectral indices 
[normalized difference spectral indices (NDSI), ratio spectral indices (RSI), different spectral indices (DSI)], and then three- 
band spectral indices were established based on the best two-band spectral index within the range of 350-2500 nm to 
reduce the noise caused by nitrogen and saturation. Then, optimal spectral indices were selected to construct models of 
LEWT monitoring in wheat. The results showed that the two-band spectral index NDSI(Ri204/ Risis) could be used for LEWT 
monitoring throughout the wheat growth season, but the model performed differently before and after anthesis. Therefore, 
further two-band spectral indices NDSIb(Ri445, R487)/ NDSIa(Ri7i4, R1395), and NDSI(Ri429, R4i6)/ were constructed for the two 
developmental phases, with NDSI(Ri429, R416) considered to be the best index. Finally, a three-band index (Ri429-R4i6~ 
Ri865)/(Ri429+R4i6+Ri865)/ whlch was superlor for monitoring LEWT and reducing the noise caused by nitrogen, was formed 
on the best two-band spectral index NDSI(Ri429, R416) by adding the 1,865 nm wavelenght as the third band. This produced 
more uniformity and stable performance compared with the two-band spectral indices in the LEWT model. The results are 
of technical significance for monitoring the water status of wheat under different nitrogen treatments in precision 
agriculture. 
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Introduction 

Real-time, non-destructive monitoring of crop water content 
based on hyperspectra is an important area of research in 
precision irrigation in agriculture [1-7]. As a widely used measure 
of crop water status, leaf-equivalent water thickness (LEWT) and 
canopy-equivalent water thickness (CEWT) not only directly 
indicate crop water content, but also provide information for leaf 
area indices. Therefore, they can visually reflect crop water 
requirements and crop growth status. CEWT has been found to 
be linearly related to the vegetation water content (VWC), with an 
value of 0.87 for corn [8]. LEWT was shown to have a better 
linear correlation with reflectance than fuel moisture content 
(FMC) in the leaves of 10 plant species [9]. In cowpeas, beans, and 
sugar beet, the R1300/R1450 leaf water index (ratio of reflectance at 
1,300 to 1,450 nm) displayed a characteristic logarithmic corre- 
lation with LEWT [10]. During the late period of wheat 
development (after anthesis), LEWT is more useful than FMC 
for assessing the water status of wheat. Several optimal water 



indices for different stages of wheat development are available 

[11]- 

Many studies in recent decades have aimed to evaluate the 
LEWT using remote sensing. Detection of plant water stress 
through remote sensing has been proposed using indices based on 
the near-infrared (NIR, 700-1,300 nm) and the middle -infrared 
(MIR, 1,300-2,500 nm) wavelengths. Hunt et al. found the 
moisture stress index (MSI) linearly correlated with the lo- 
gioLEWT of Quercus agrifolia (sclerophyllous leaves) and Liquidambar 
styracijlua (hardwood deciduous tree leaves), in which the regression 
equations were different [12]. Ceccato et al. reported that 
shortwave infrared (SWIR, 1,400-2,500 nm) was sensitive to 
LEWT, but could not be used alone to determine LEWT because 
two other leaf parameters (internal structure and dry matter) also 
influence leaf reflectance in the SWIR [13]. A combination of 
SWIR and NIR was necessary to determine LEWT. Gao 
proposed a new vegetation index, the normalized difference water 
index (NDWI), for the remote sensing of CEWT from space [14]. 
This index was constructed based on two narrow bands centered 
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Table 1. Basic information about different field experiments. 



Treatment (water 
treatments (W): %; 

Year and nitrogen rates Spectrum and Soil (yellow Environmental conditions 

Exp. Number Variety (N): kg/hm^)) sample data brown soil) during the growing seasons 



Exp. 1 (Calibration 
data set) 



Exp. 2 (Validation 
data set) 



2010-2011 
Yangmai 18 



2011-2012 
Yangmai 18 



4 W: W1 (9.5-10.5), 
W2 (15.5-16.5), W3 
(21.5-22.5), W4 
(29.5-30.5). 2 N: 
N1 (150), N2 (300) 

3 W: W1 (13.5-14.5), 
W2 (21.5-22.5), W3 
(29.5-30.5). 3 N: N1 
(90), N2 (180), N3 (270) 



Jointing (3.25), 
Booting (4.1), 
Heading (4.12), 
Anthesis (4.20), 
Filling (5.6) 

Jointing (3.24), 
Booting (4.3), 
Heading (4.11), 
Anthesis (4.16), 
Filling (4.22) 



Organic matter: 15.5 g kg \ 
Total N: 1.1 g kg \ Available 
P: 50.8 mg kg \ Available K: 
89.6 mg kg \ 

Organic matter: 14.8 g kg \ 
Total N: 1.1 g kg \ Available 
P: 50.4 mg kg \ Available K: 
88.9 mg kg \ 



Mean temperature: 10.84°C Maximum 
temperature:1 5.75Minimum 
temperature:5.92Mean Diurnal temperature: 
9.83°C Sunshine hours: 1302.5 
hPrecipitation:1 55.14 mm 

Mean temperature: 10.98°C Maximum 
temperature:1 4.92Minimum 
temperature:7.04Mean Diurnal temperature: 
7.88°C Sunshine hours: 999.003 
hPrecipitation:329.51 mm 
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near 860 nm and 1,240 nm and was used successfully to detect 
GEWT and LEWT in cotton and trees [15-18]. Based on the 
difference in reflectance between 945 nm and 975 nm and using 
Beer's law, Liu et al. calculated the radiation-equivalent water 
thickness of leaves (RLEWT) [19]. The authors demonstrated that 
RLEWT was significantly correlated with LEWT. Zarco et al. 
estimated LEWT from canopy-level reflectance with the simple 
ratio water index (SRWI), which had a strong correlation with 
LEWT [20]. Other researchers have proposed the normalized 
difference infrared index (NDII) [8,21] and the water index (WI) 
[22] to estimate LEWT and CEWT. 

In addition to these two-band spectral indices, three-band 
spectral indices have been proposed to evaluate other growth 
parameters of plants. Schneider et al. [23] and Stow et al. [24] 
found the visible atmospherically resistant index (VARI) to be 
minimally sensitive to atmospheric effects and strongly related to 
live fuel moisture (LFM). Li et al. [25] and Jie et al. [26] found 
VARI- 7 00 to be significantly correlated with yield at the whole 
development stages in cotton. Wang et al. constructed three-band 
vegetation indices, (Ra.i-R^2+2 xR^3)/(R^i-R:>,2~2 xR^^s) and 
(Rp^i— R;,2~R?^3)/(R?^i+R?^2+R?^3)5 to reduce the saturation ob- 
served in two-band vegetation indices [27]. They demonstrated 
that the models for leaf nitrogen content (LNC) using (R924~ 
R703+2xR423)/(R924+R703~2xR423) wcrc Stable and accurate 
and more effective than other published vegetation indices. 

Since the value of R445 is constant until the total chlorophyll 
content drops below 0.04 mol m ^, Sims et al. added R445 to 
ND705 and SR705 and developed new three-band spectral indices. 
The modified ND705 and SR705 (mND705, mSR705) were used to 
predict leaf pigment content, and both were insensitive to species 
and leaf structure variation [28]. Tian et al. developed a blue 
nitrogen index (R434/(R495+R4oi)) and estimated the canopy leaf 
nitrogen concentration of rice [29] . 

Water and nitrogen are the main limiting factors in plant 
growth [30-32], and they interact in a complex manner. Previous 
studies have shown that water and nitrogen are indirectly related 
(within plants) through chlorophyll and cellulose [33]. It may be 
effective to construct spectral indices that include a waveband 
sensitive to chlorophyll or cellulose, which might indirectly 
eliminate the impact of nitrogen on LEWT monitoring with good 
performance. 

In this study, two experiments were conducted on winter wheat 
with different water and nitrogen treatments in two consecutive 
growing seasons. The objectives of the study were: (1) to determine 
hyperspectral bands that were sensitive to LEWT but insensitive to 



nitrogen in wheat, (2) to develop new spectral indices for 
monitoring LEWT in wheat, and (3) to quantify the relationships 
between LEWT and the new spectral indices to reliably estimate 
LEWT. These results may provide a technical approach to 
effectively monitoring of plant water status while minimizing the 
noise from nitrogen in precision wheat management. 

Materials and Methods 

2.1 Design of field experiments 

Two experiments were conducted at the Pailou Experiment 
Station at Nanjing Agricultural University, China (118°15'E, 
32°rN). These experiments involved different water (W) and 
nitrogen (N) treatments. The wheat cultivar Yangmai 18 was 
grown in two consecutive seasons from November 2010 to June 

2011 (a low rainfall season) and from November 2011 to June 

2012 (a high rainfall season). Seeds were sown on 1 November 
2010 and 3 November 201 1 at a density of 180 plants per m^ with 
a plot size of 10 m^ (2.5 mx4 m). Each experiment employed a 
randomized complete block design with three replications. The 
data from Experiment 1 were used to derive the monitoring 
models, while the data from Experiment 2 were used to evaluate 
the models. More details about the W and N treatments, sampling 
procedures, and environmental conditions are given in Table 1. 

There were 24 plots in 2010-201 land 27 in 2011-2012. Each 
plot was constructed from cement with identical length (3 m), 
width (3 m) and depth (1 m). A transparent plastic sheet was 
placed over the plots to a height of 3 m above the ground to 
prevent natural rainfall from flowing into the plot and disturbing 
the experimental water treatments. From greening until harvest, 
the runoff from the plastic covers was discharged into a cement 
drainage channel outside of the plots. From early jointing, the 
volumetric soil water content was measured with a TRIME-PICO 
TDR (TRIME-PICO, IMKO, Germany) portable soil moisture 
speed measuring device in each plot using the five-point method at 
4:00 p.m. local time every day. If necessary, water was added 
immediately to the soil to maintain the designed water content. 
The date for obtaining spectrum and plant samples was jointing, 
booting, anthesis, filling. 

2.2. Data measurements 

2.2.1. Measurement of leaf hyperspectral 
reflectance. Leaf spectral measurements were taken using an 
accessory of the ASD Field Spec Pro spectrometer (Analytical 
Spectral Devices, Boulder, CO, USA). The accessory comprises a 
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handheld leaf folder spectral detector with its own light source, 
which is designed to reduce the effect of time of day, weather, 
atmospheric vapor pressure, or soil background on the readings. 
The ASD spectrometer is operated in the 350-2,500 nm spectral 
region, with a sampling interval of 1.4 nm and spectral resolution 
of 3 nm between 350 and 1,050 nm, and a sampling interval of 
2 nm and spectral resolution of 10 nm between 1,050 and 
2,500 nm. 

Ten wheat plants were randomly sampled from each plot. The 
top four leaves on each plant were identified (numbered LI to L4 
starting from the top), and their spectral reflectance was measured. 
The spectral reflectance averaged over 10 plants at each leaf 
position was taken as the reflectance for that leaf position. Before 
the measurement of leaf reflectance in each plot, a standard 
whiteboard (Labsphere, North Sutton, NH, USA) was used to 
calibrate the spectral reflectance of the leaves. 

2.2.2. Determination of leaf equivalent water 
thickness. After the leaves were measured for spectral reflec- 
tance, they were stored in pre-weighed valve bags in ajar of liquid 
nitrogen. The leaves were transferred to the laboratory for the 
determination of leaf area using a portable leaf area meter (LAI- 
3000, Licor, NE, USA) and their fresh and dry weights were 
obtained. LEWT was calculated as follows: 

LEWT = (WF-WD)/(DW * A) * 10000(|im), 

where Wf is the leaf fresh weight (g), Wd is the leaf dry weight (g), 
Dw is the water density value (g/ cm^), and A is the leaf area (cm^). 

2.2.3 Determination of leaf nitrogen content. For each 
sample, the leaf dry weight was determined by oven-drying the 
leaves at 80°C to a constant weight. The LNC was determined on 
a dry weight basis (g 100 g ^) using the micro-Kjeldahl method. 

2.3. Data analysis 

The LEWT values and the corresponding leaf spectral 
reflectance for each sampling date were analyzed with MATLAB 
8.2 and Excel 2011 software. The two-band spectral indices 
considered included the normalized difference spectral index 
(NDSI, equation (1)), the ratio spectral index (RSI, equation (2)), 
and the difference spectral index (DSI, equation (3)). To determine 
the wavelengths in the two-band spectral indices, all combinations 
of wavelengths within the range of 350—2,500 nm at 1 nm 
intervals were evaluated according to the criteria presented below: 

NDSI = (Rai -Ra2)/(Rai + Ra2) ( 1 ) 



RSI = R,i/Ra2 (2) 



DSI = R,i-R,2, (3) 

where Rp^i is spectral reflectance at the wavelength of ^i. 

For the three-band spectral indices, the best wavelengths were 
combined with all possible third wavelengths within the range of 
350—2,500 nm at 1 nm intervals. The principle of retaining the 
best wavebands of the three-band indices has been applied in 
previous studies [26-27]. Three classical three-band spectral 
indices based on the NDSI were considered and constructed 
(equations (4) to (6)). 

(Rai-Ra2)/(Rai + Ra2-Ra3) (4) 



(R^1-R;12-Ra3)/(Ra1 + Ra2 + Ra3) (5) 



(Rai-Ra2 + 2R^3)/(Rai + Ra2-2Ra3) (6) 

The criteria for evaluating the best wavelengths were the 
coefficient of determination (R^), standard error (SE), and relative 
root mean square error (RRMSE) [25]. 

The RRMSE was calculated as follows: 



/l " 100 
RRMSE = J-^Y. - ^')' ^ 

where Pi and 0^ are the predicted and observed values, 
respectively, 0/ is the observed mean value, and n is the number 
of samples. 

Results 

3.1 Variation in LEWT 

Fig. lA shows how LEWT varied under the same nitrogen level 
and different water treatments in Exp. 1. There was a greater 
difference at anthesis and filling than at booting, jointing, or 
heading due to water treatment, with LEWT values in the order of 
W4 > W3 > W2 > Wl, which is normal for irrigation scenarios. 
During the late stages of plant development, evaporation increased 
with the air temperature, which resulted in greater differences in 
LEWT among water treatments. Drought stress (Wl, W2) was 
found to accelerate plant development, and LEWT value was 
higher than in W3 (normal amounts of water), while the W4 
treatment was found to enhance the acceleration of plant 
development. 

Fig. IB displays LEWT for the W3 treatment for leaf LI at 
various nitrogen levels. LEWT was greater at N2 than at Nl, 
which was also the case for the Wl, W2, and W4 treatments. Due 
to the synergistic effects of water and nitrogen on the growth 
processes of wheat, high-level nitrogen treatment promoted water 
uptake, resulting in a greater LEWT compared to low nitrogen 
levels [32-33]. Similar results were obtained for L2, L3, and L4 in 
Exp. 2. 

3.2 Variations in leaf hyperspectral reflectance 

Differences in LEWT and LNC values significantly affected the 
leaf hyperspectral reflectance of wheat. Fig. 2A shows leaf 
reflectance at different LEWT levels for N2 during the booting 
stage in Exp. 1. Leaf reflectance in the visible range was not 
significantly aflected by LEWT. However, in the NIR range, 
especially at the central wavelengths of 900 nm, 1,200 nm, 
1,400 nm, 1,450 nm, and 1,930 nm, leaf reflectance significantly 
decreased as LEWT increased. This result suggests that nitrogen 
levels influenced the variations in leaf hyperspectral reflectance. 
Thus the reflectance spectrum used for monitoring crop water 
may contain some information due to the presence of nitrogen, 
which should be eliminated or reduced. 

Fig. 2B shows the reflectance at various LNC levels for the W4 
treatment. N levels and W treatments markedly influenced the 
characteristics of spectral reflectance, with different spectral 
responses in the various waveband regions. With increasing N 
levels, the reflectance decreased in the visible bands and increased 
in the near infrared, with obvious differences among the four N 
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Figure 1. Changes of LEWT under different water treatments at N2 (A) and varied nitrogen treatments at W3 (B) in Exp. 1. 
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levels. This result implies that these spectral regions are relatively 
sensitive to the growth status of wheat at different N levels. In the 
visible range, differences in reflectance appeared mainly at 
wavelengths of 450 nm and 700 nm, the central wavelengths of 
pigment absorption. The reflectance increased as LNC increased 
in the near infrared, particularly at 970 nm and 1,200 nm, the 
central wavelengths of water absorption. 

3.3 Correlation of LEWT and LNC with the original and 
first derivative hyperspectra in wheat leaves 

The correlation of LEWT and LNC with the original spectrum 
and its first derivative spectrum was determined based on the 
observed data from Exp.l. Fig. 3 shows that LEWT and the 
original spectral reflectance were positively correlated in the visible 
region and negatively correlated in the NIR, with peaks of the 
absolute value of the correlation at 970 nm, 1,200 nm, 1,395 nm, 
1,450 nm, 1,850 nm, and 2,200 nm. The strongest correlation 
was at 1,395 nm (R1395 = —0.7576). This result indicates that 
LEWT markedly influenced the reflectance of wheat leaves, with 



different spectral responses in the various waveband regions. In 
general, the correlation between LNC and the original spectral 
reflectance generally had the opposite sign to the correlation of 
LEWT. 

Fig. 3 A shows that the largest 10% of |R|lewt ^^d |R|lnc 
overlapped at 1,365-1,424 nm and at 1,839-1,869 nm in the NIR 
region. The largest 20% of |Rlewt| and |Rlnc| overlapped in 
two regions of the visible bands, 526-566 nm and 701-723 nm 
(Fig. 3A). These results indicate that the wavelengths sensitive to 
both LEWT and LNC were mainly in the visible and NIR range. 
However, the largest 20% of |Rlewt| did not overlap anywhere, 
with the smallest values of | Rlnc I •> suggesting that no wavelength 
was sensitive to LEWT and insensitive to LNC in the same region. 

The first derivatives of the original spectrum and LEWT were 
also strongly correlated in the NIR range, especially at the central 
wavelengths of 1,890 nm, 1,300 nm, and 1,400 nm (Fig. 3B). The 
largest 10% of |Rlewt| overlapped with | Rlnc I 770 nm 
(Fig. 3B). In the original spectrum, there was no overlap between 
the largest 20% of |Rlewt| and the smallest values of | Rlnc I- 
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Figure 3. The correlation of LEWT and LNC to original spectral reflectance (A) and its first derivative (B) in Exp. 1. 
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No wavelength in the first derivative was sensitive to LEWT and 
insensitive to LNC simutaneously. 

3.4 The relationships of spectral indices to LEWT 

3.4.1. Relating two-band spectral indices to LEWT during 
the entire growth period. Contour maps of values for the 
relationship of LEWT during the full growth period with the two 
wavelengths in the two-band indices were constructed (not shown). 
The results for NDSI and RSI were similar, while R^ values were 
lower for DSL The top 10% of R^ values for NDSI, usins the 
calibration data, were mainly in the ranges of XI = 1,000- 
1,320 nm and X2 = 1,300-1,350 nm (Fig. 4A). Using the valida- 
tion data, the top 10% of R^ values were in the ranges of 
X\ = 1,205-1,315 nm and X2 = 1,287-1,337 nm (Fig. 4B). Finally, 
NDSI (Ri2045 Risis) was selected based on the overlapping region. 
This index had an R^ value of 0.708 and an SE of 15.5053 |im 
(Fig. 5A). Fig. 5B shows the relationship between the predicted and 
observed LEWT. There are clear differences in the fit to the data 
before and after anthesis, which may be due to metabolism 



(especially the nitrogen content) and changes in position of the 
plants. 

Further statistical analysis was conducted to determine whether 
the two linear relationships could be grouped. Thereby, a linear 
parallel curve analysis with grouped data was used to determine 
whether the critical curves for LEWT differed before and after 
anthesis [34] . The results showed that the F value for the LEWT 
before and after anthesis was significantly different from that of the 
residual. The F value for NDSI(Ri204, Ri3i8)xLEWT before 
anthesis and NDSI(Ri2045 Ri3i8)xLEWT after anthesis also 
differed significantly from that of the residual (Table 2). 

3.4.2. Relating NDSI (Rgi, R^z) to LEWT before and after 
anthesis. Based on the above results, the developmental stages 
were considered separately in further studies of the relationship 
between LEWT and the NDSI (R^i, Re2). The optimal 
wavelengths for the period before anthesis were 1,445 and 487, 
giving the spectral index NDSIb(Ri445, R487) (NDSI before 
anthesis). This index had an R^ value of 0.8137 and an SE of 
10.6503 |im. The optimal wavelengths for the period after anthesis 
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Figure 4. Contour map of coefficients of determination (R^) for linear relationship between NDSI and LEWT (A: calibration dataset, 
B: validation dataset). 

doi:1 0.1 371 /journal. pone.0096352.g004 



PLOS ONE I www.plosone.org 



5 



June 2014 | Volume 9 | Issue 6 | e96352 



Exploring Novel Bands and Key Index for LEWT in Wheat 




Figure 5. Quantitative relationships between NDSI (R1204/ Risis) and LEWT (A); the 1:1 relationship between the predicted and 
observed LEWT based on NDSI (R1204/ Risis) (B). 
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were 1,714 and 1,395, giving the spectral index NDSIa(Ri7i4, 
R1395) (NDSI after anthesis). For this index, the value was 
0.8622 and the SE was 12.2442 jim (Fig. 6). The independent 
dataset from Exp. 2 was used to test the models for LEWT before 
and after anthesis (Fig. 7). The results showed that NDSIb and 
NDSIa predicted LEWT more accurately than did the spectral 
index NDSI (R1204? Risis) for the entire growth period and the 
spectral index NDSI(Ri429, R416) of the common sensitive areas. 
In the optimal spectral indices NDSIb(Ri445, R487) and 
NDSIa(Ri7i4, R1395), 1,445 nm and 1,395 nm were the most 
sensitive wavelengths to LEWT, as explained above. The bands at 
487 nm and 1,714 nm were sensitive to chlorophyll and lignin. 

The indices NDSIb(Ri445, R487) and NDSIa(Ri7i4, R1395) 
contain four wavelengths, which will increase the cost of 
manufacturing-self This motivated us to search for common 
two-band wavelengths before and after anthesis. The central 
waveband ranges of NDSI before and after anthesis differed 
significantly, and the top 10% of the optimal values of R^ had no 
common area (Figures were omitted). However, the top 15% of R^ 
values before (Fig. 8A) and after (Fig. 8B) anthesis had two 
common areas, with the first X\ in the range of 405-420 nm and 
X2 in the range of 1,418-1,477 nm and the second XI in the range 
of 408-443 nm and X2 in the range of 1,876-1,886 nm (Fig. 8C). 
Based on the R^, SE, and RRMSE of the calibration and 
validation data, the spectral index NDSI(Ri429, R416) was selected 
for the periods both before and after anthesis. The calibration of 
LEWT both before and after anthesis had an R^ value of 0.6776 



and an SE of 16.2914, and the prediction using the validation data 
had an R^ value of 0.3967 and an RRMSE of 0.1243 (Fig. 9, 
Table 3). 

3.4.3. Relating three-band spectral indices based on 
NDSI(Ri429, R416) to LEWT. The third band {X3) can help to 
reduce both the effect of nitrogen and the effect of saturation. The 
best three-band index, using the best wavelengths for NDSI, was 
(Ri429-R4i6-Ri865)/(Ri429+R4i6+Ri865)- The third band (at 
1,865 nm) located in the SWIR is close to the sensitive band of 
cellulose (1,880 nm) (Table 3). However, the performance of the 
model for the three-band index for NDSIb(Ri445, R487) and 
NDSIa(Ri7i4, R1395). The cost of developing the equipment in the 
future will decrease as its use on a large scale becomes more 
feasible due to the need for fewer wavelengths. Moreover, the 
three-band index model displayed a marked improvement 
compared with the model based on NDSI(Ri429, R416) (Fig. 9, 10). 

3.4.4. Comparison of the performance of the new and 
previous models. Five popular existing spectral indices that are 
used to predict crop water content were compared with the 
spectral indices developed in the current study. As shown in 
Table 3, of the five indices, WI, MSI, and NDII yielded the best 
predictions of LEWT, and WI performed especially well. 
However, the two spectral indices developed in this study, NDSIb 
and NDSIa, performed considerably better than all of the existing 
indices. This result may be due to the fact that the existing spectral 
indices are based on the relative water content of leaves and whole 
plants of various crops, whereas the spectral indices developed in 
the current study were constructed based on the LEWT in wheat. 



Table 2. Simple grouping linear analysis of the model between NDSI and LEWT in wheat before and after anthesis. 





Developmental stage 




NDSI (Ri204/Ri3i8) xdevelopmental stage 


Residual 




df MS 


F 


df MS F 


df MS 


F 



1 7478.7 37.34^ 1 4202.6 20.98^ 267 200.3 6.76 



Notes: 

df: degrees of freedom. 
MS: mean square. 
^: P<0.05. 

doi:1 0.1 371/journal.pone.0096352.t002 
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Figure 6. Contour map of the linear relationship with top 1 5% between NDSI and LEWT before (A) and after (B) anthesis, and the 
common area between A and B (C). 

doi:1 0.1 371 /journal. pone.0096352.g006 
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Figure 7. Quantitative relationships between NDSIe (R1445, R487) and LEWT before anthesis (A); quantitative relationships between 
NDSIt (Ri7i4, R1395) to LEWT after anthesis (B). 

doi:1 0.1 371 /journal. pone.0096352.g007 
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Figure 8. The 1:1 relationship between the predicted and observed LEWT based on NDSIb (R1445, R487) for wheat before anthesis 
(A); the 1:1 relationship between the predicted and observed LEWT based on NDSIa (R1714, R1395) for wheat after anthesis (B). 

doi:10.1371/journal.pone.0096352.g008 
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Figure 9. Quantitative relationships between (Ri 

429~R4i6)/(Ri429+R4i6) 3nd LEWT (A); the 1:1 relationship between the predicted 
and observed LEWT on (Ri429 — R4i6)/(Ri429"l"R4i6) wheat (B), 

doi:1 0.1 371 /journal. pone.0096352.g009 



Discussion 

4.1. Response of leaf spectral reflectance to water and 
nitrogen treatments in wheat 

Previous studies of the response of leaf spectra to different water 
conditions and nitrogen levels have shown that the sensitivity of 
leaf spectral reflectance to water content peaks at the central 
wavelengths of 970 nm, 1,200 nm, 1,400 nm, 1,450 nm, 
1,730 nm, 1,930 nm, 2,100 nm, and 2,500 nm [35]. Within the 
range of 1,300-2,500 nm, spectral reflectance is affected by water 
in the plants, which directly absorbs radiation. Within the range of 
400-1,300 nm, reflectance is influenced by changes in the internal 
structure of the leaf, which is a result of their different water 
content ratios [7,36]. The influence of nitrogen on spectral 
reflectance occurs mainly in the range of the visible and NIR 
wavebands; and the wavebands most sensitive to nitrogen are in 
the range of 530-560 nm [37-38]. In the present study, leaf 
hyperspectral reflectance at different LEWT levels in the N2 
treatment did not change significantly in the visible band range. 
However, reflectance significantly decreased as LEWT increased 
within the NIR range, especially at 900 nm, 1,200 nm, 1,400 nm, 
1,450 nm, 1,730 nm, and 1,930 nm. Leaf hyperspectral reflec- 
tance increased with increasing LNC within the NIR range, and it 
decreased considerably at the central wavelengths of water 
absorption. Therefore, to monitor plant water content, some 
sensitive wavelength and spectral indices that are strongly 
correlated with LEWT, but less well correlated with LNC, should 
be extracted with consideration of the fact that the spectrum 
reflectance is affected by both water and nitrogen. 

4.2. Two-band spectral indices sensitive to LEWT in wheat 

The systematic analysis of two-band spectral indices sensitive to 
LEWT was conducted for the entire growth. NDSI is more closely 
related to LEWT than RSI or DSL NDSI(Ri204, Risis) was the 
optimal spectral index for predicting LEWT in wheat during the 
whole growth period. As earlier analyses have indicated, 
differences in LNC before and after anthesis might have an 
impact on the accuracy of the LEWT model throughout plant 
growth. Therefore, separate indices were proposed for the periods 
before and after anthesis, namely, NDSIb(Ri445, R487) and 



NDSIa(Ri7i4, R1395), respectively. The wavelengths 1,445 nm 
and 1,395 nm were the most sensitive to LEWT, while the 
wavelengths 487 nm and 1,714 nm were sensitive to chlorophyll 
and lignin, respectively. This indicated that LEWT models before 
and after anthesis were established separately, while the impact of 
nitrogen was eliminated by the use of wavebands that were 
sensitive to chlorophyll or lignin, as reported previous [33]. 

4.3. Three-band spectral indices sensitive to LEWT in 
wheat 

Three-band indices contain less noise due to chlorophyll or 
lignin than two-band spectral indices [23-24,28]. They also 
display less saturation than two-band spectral indices [27]. 
Therefore, to reduce the noise caused by nitrogen before and 
after anthesis in the LEWT model, the three-band spectral index 
(Ri429-R4i6-Ri865)/(Ri429+R4i6+Ri865) was developed bascd 
on the best overall two-band spectral index NDSI(Ri429, R4i6)- 
Although the model based on the three-band index was a poorer 
predictor than those based on NDSI (R12045 Risis) or on 
NDSIb(Ri445, R487) and NDSIa(Ri7i4, R1395), the calibration 
and validation on the resulting scatter plot was centralized, and the 
LEWT model was more uniform and stable with a lower RRMSE. 
The results also showed that the three-band spectral indices 
(R 1 429 - R4 1 6 - R 1 865) / (R 1 429+R4 1 6+R 1 865) clearly reduced the 
noise due to nitrogen in the LEWT model throughout the whole 
period of wheat growth, with an R^ value of 0.7255and RRMSE 
of 0.0689 for LEWT validation, and an R^ value of 0.0849 and an 
RRMSE of 39.7950 for LNC validation. 

Conclusions 

Previous studies have shown that the transportation of amino 
acids from leaves to grains after anthesis leads to physiological and 
biochemical changes in the organizational structures of leaves, 
which affect LEWT monitoring based on leaf hyperspectral 
reflectance. In this study, we demonstrated that different water 
and nitrogen treatments affected the variation in LEWT and the 
leaf hyperspectral reflectance of wheat within the 350-2,500 nm 
range. Furthermore, the top 10% of the maximum |Rlewt| and 
I Rlnc I values were found to share common wavelength ranges. 
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Figure 10. Quantitative relationships between (Ri429-R4i6-Ri865)/(Ri429+R4i6+Ri865) between LEWT in wheat (A); the 1:1 
relationship between the predicted and observed LEWT on (Ri429~R4i6~Ri865)/(Ri429+R4i6+Ri865) (B). 

doi:10.1371/journal.pone.0096352.g010 



Based on this study and previous reports, when monitoring 
LEWT, the noise from the LNC should be considered. The model 
based on a three-band index (Ri429~R4i6~Ri865)/(Ri429+R4i6+ 
R1865) described in this paper performed well for monitoring 
LEWT, with a higher predictability and stability for water content 
and lower noise levels due to nitrogen under various water and 
nitrogen treatments. 
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